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Abstract: We study the moduli space oi 4d M = 1 supersymmetric QCD in the Veneziano 
hmit using Hilbert series. In this hmit, the numbers of colours and flavours are taken to be 
large with their ratio fixed. It is shown that the Hilbert series, which is a partition function 
of an ensemble of gauge invariant quantities parametrising the moduli space, can also be 
realised as a partition function of a system of interacting Coulomb gas in two dimensions. 
In the electrostatic equilibrium, exact and asymptotic analyses reveal that such a system 
exhibits two possible phases. Physical quantities, such as charge densities, free energies, and 
Hilbert series, associated with each phase, are computed explicitly and discussed in detail. 
We then demonstrate the existence of the third order phase transition in this system. 



Contents 



1 Introduction and summary 2 

2 The Hilbert series of SQCD and the Coulomb gas 5 

2.1 The continuum limit 6 

3 Exact results in the continuum limit 8 

3.1 The charge density of the Coulomb gas 8 

3.2 The integration limits in the Coulomb gas 10 

3.3 The free energy 12 

4 Series expansions in the Veneziano limit 14 

4.1 The gapless phase 15 

4.2 The gapped phase 17 

5 Leading order approach in the Veneziano limit 20 

5.1 Leading order solution 21 

5.2 Linear perturbation to the free energy 24 

5.3 Discussion 26 

6 Phase transition 26 

6.1 The gap width near the transition point 27 

6.2 Free energy at the phase transition 27 

A Integral identities 29 

A.l Integrals for the charge density 29 

A. 2 Integrals for the free energy 29 

B Method of conformal mappings and mirror charges 32 

B. l Solving the complex potential 32 
B.2 The free energy 33 
B.3 The solution for the gapless phase 34 
B.4 The solution for the gapped phase 34 
B.5 Charge densities 35 



- 1 - 



1 Introduction and summary 



Hilbert series of a supersymmetric gauge theory is a partition function of the system of 
gauge invariant quantities parametrising the moduH space of vacua. It can be viewed as 
a trace over the space of gauge invariant objects parametrising the vacua of the theory, 
counting with respect to a certain U{1) global symmetry [1-9]. In the context of 4d A/" = 1 
supersymmetric QCD (SQCD), these gauge invariant objects are combinations of quark and 
antiquark superfields, known as mesons and baryons [6, 7]. It is convenient to take the U{1) 
global symmetry to be proportional to the i?-symmetry such that the quarks and antiquarks 
carry a unit charge. 

In computing the partition function, we need to introduce a bookkeeping variable of the 
C/(l) global charge; this is known as the fugacity denoted by t, such that < t < 1. The 
Hilbert series can be computed in two steps: 

1. Compute the generating function of symmetric functions of the quark and antiquark 
superfields. This is a rational function of t and Cartan variables of the gauge group G, 
denoted by for a = 1, 2, . . . , rank G. 

2. Integrate this rational function over the Haar measure of the gauge group G. This step 
restricts the generating function in Step 1 to count the gauge invariant combinations of 
quarks and antiquarks, namely mesons and baryons. 

Since the Haar measure involves integrals along the contours \za\ = 1 for a = 1,2,..., rank G, 
the Hilbert series can be written in terms of nested multi-contour integrals over a certain 
rational function. These expressions are explicitly written in (4.9) of [6], (3.14) of [7], (3.6) 
of [10], (2.8), (3.3), (4.4) of [11], and (5), (8) of [12]. 

In [10], [11], and [12], it was shown that these multi-contour integrals can be recast in 
terms of determinants of Toeplitz matrices for unitary and special unitary gauge groups, and 
determinants of Hankel matrices for the special orthogonal and symplectic gauge groups. 
Writing these integrals in terms of determinants has several advantages: it allows for exact 
expressions for the Hilbert series for arbitrary number of colours and flavours, as well as 
various asymptotic results when the number of colours and flavours are large. 

One of the asymptotic limits considered in this paper is when the number of colours Nc 
and the number of flavours Nf tend to infinity, with their ratio kept fixed and finite. This 
is known as the Veneziano limit [13]. Historically, this limit has received a lot of interest in 
QCD due to the interesting topological nature of Feynman diagrams. The main goal of this 
paper is to study the vacuum moduli space of SQCD via analyzing the Hilbert series in this 
asymptotic limit. 

Another motivation to study the Veneziano limit of SQCD is that given the ratio A'^^ /Nc 
in a certain range, the theory may exhibit a number of physically interesting aspects. Let us 
draw some examples from [14-16] as follows: 

• For SU(Nc), SO{Nc), and Sp{Nc) gauge groups, with respectively Nf>Nc + 2, Nf > 
Nc + 3, and Nj > A'c + 3, each of these theories not only possesses a degenerate vacua at 
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the quantum level, but its classical moduli space also receives no quantum corrections. 
Hence, the Hilbert series computed for the classical moduli space as described above 
still provides a valid description of the quantum moduli space. 

• In these ranges of Nf and Nc, the theory possesses a dual description. This phenomenon 
is known as the Seiberg duality [17]. The Hilbert series will therefore also describe the 
quantum moduli space of the dual theory. 

• If Nf/Nc is fixed in a certain interval, known as the conformal window, the theory has 
a non-trivial infrared fixed point. For SU{Nc), Sp{Nc), and SO{Nc) gauge groups, 
the conformal windows are §7Vc < Nf < 3Nc, |(iVc + 1) < iV/ < 3(iVc + 1), and 
|(A^c ~ 2) < Nf < 3{Nc — 2), respectively. Within this interval, the Hilbert series 
describes the vacua of the theory at the conformal fixed point. 

The moduli space of U{Nc) and SU {Nc) SQCD with Nf flavours in this asymptotic limit 
was studied in detail in [10, 12]. In this paper, we focus on SQCD with SO{Nc) and Sp{Nc) 
gauge groups. 

Let us now present an overview of the approach we use in this paper. Regarding the 
Hilbert series as a partition function of the system of gauge invariant objects, we can read 
off the Hamiltonian of this system from the multi-integral form of the Hilbert series. Such 
a Hamiltonian in the Veneziano limit can be interpreted as the Hamiltonian of a system of 
Coulomb gas or log- gas in two dimensions. The charged particles are constrained to reside 
in a one-dimensional interval inside the two-dimensional space. The charges in this interval 
interact among themselves and are also subject to the interactions due to external fields; 
the latter is determined by the number of ratio Nf/Nc and the gauge group. By solving 
the equilibrium condition for this system, we can compute the charge density within the 
aforementioned one-dimensional interval. This quantity exhibits two phases of the system, 
dubbed the gapless and the gapped phases, with the order parameter being the size of the gap. 
Which of these two phases the system belongs depends on the ratio Nf /Nc and on the value 
of the fugacity t in the Hilbert series. The free energy of the system can be computed in each 
phase, from which the Hilbert series can be deduced. By examining the behavior of the free 
energy at the transition point between these two phases, we demonstrate the possibility of 
the third order phase transition, reminiscient of [18, 19]. This phase transition was uncovered 
for the U{Nc) and SU{Nc) gauge groups in [12]. 

Outline and key results 

This paper is divided into four main parts. The first part, containing Section 2, develops 
the connection between Hilbert series and the partition function of the Coulomb gas in two 
dimensions, which was established for the U{Nc) and SU{Nc) gauge groups in [12]. The 

^It was shown in [14-16] that, for SU{Nc) SQCD with Nf < - 1, SO{Nc) SQCD with Nf < Nc - 5, 
and Sp{Nc) SQCD with Nf < Nc, the moduh spaces are totally lifted by dynamical generated superpotential 
and hence there are no supersymmetric vacua. 
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basic connection goes back to the seminal works of Wigner and Dyson [20-25]. While the 
electrostatic analogues have a long history in string theory [26, 27] (see also [28, 29]), they 
have received considerable attention only relatively recently (see e.g., [30-44]). By exploiting 
the Coulomb gas picture we coarse-grain the Hamiltonian of the Coulomb gas system in the 
Veneziano limit and determine the electrostatic equilibrium condition. 

We then proceed to the second part of the paper, consisting of Sections 3 and 4. In 
Section 3, we generalise the method of [45] to find exact solutions to this equilibrium condition, 
from which expressions for the charge density and the free energy are derived. These exact 
expressions involve certain variables determined by degree six polynomial equations. Thus, 
in order to make these exact expressions explicit in the Veneziano limit, we need to solve 
such equations perturbatively in 1/Nc (with the ratio Nf/Nc kept fixed) and this is the main 
goal of Section 4. These perturbative solutions suggest that the Coulomb gas system at the 
equilibrium exhibits two phases. The charge density, the free energy, and the Hilbert series 
for each phase are computed. The agreement between the Hilbert series computed in this 
paper and known results (or limits of known results) presented in [7, 11] justifies the validity 
of our perturbative solutions. 

In the third part of the paper, namely Section 5, we present another approach to derive 
the key results for the Hilbert series. We apply directly a perturbative analysis to the Coulomb 
gas system in the Veneziano limit. We solve the electrostatic equilibrium condition at leading 
order - the next to leading order correction to the free energy, and consequently to the Hilbert 
series, then immediately follow by the standard perturbation theory approach. This is enough 
to reproduce the main results of Section 4. 

In the last part of the paper, we examine the behaviour of the Coulomb gas system near 
the transition point between the two phases. The order parameter and the free energy are 
investigated. We demonstrate the situation in which the system exhibits the third order phase 
transition. 

We summarise the key points of this paper below. 

• The two phases of the system are discussed and summarised above (4.1) and below 
(5.5). The equilibrium configurations of the Coulomb gas for each phase are solved 
explicitly in Sections 4.1, 4.2, and 5. 

• The Hilbert series computed in the gapless phase are given by (4.13) and those computed 
in the gapped phase (in a certain limit) are given by (4.29)-(4.31). 

• The existence of the third order phase transition is shown in (6.8). 
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2 The Hilbert series of SQCD and the Coulomb gas 



The Hilbert series for SQCD with classical gauge groups SO{2n + 1), S0{2n), and Sp{n) 
with Nf flavours can be written as [11] 



9Nf,G{t) = < 



_2 



(l-t)"/(2t) 

2(n-l) +n 



j^Vcin, Nf) for G = SO{2n + I) = B„ 



for G = SO{2n) = Dn 
for G = Sp{n) = Cn , 



where Vein, Nf) is defined by 

VG{n,Nf) = dxi--- J dxn n (^'^ 

^ l<a<b<n 



n 



Wg{Xc) 



-Xc + 



1+12 
2t 



(27r)"n! 



Y[ {Xa - Xb 

(-1-1)" l<a<b<n c=l 



where 



T :-- 



2t 



with T > 1 



and wg{x) is known as the Jacobi weight defined by 

wg{x) = (l-x)"(l + x)^ , 

with 



(a, 13) 



(1/2,-1/2) for G = SOi2n + I) = B„ 
(-1/2,-1/2) ior G = SO{2n) = Dn 
(1/2,1/2) for G = Sp{n) = Cn ■ 



(2.1) 



n«^G(xc)e-^/i°s(^-^'^)dx, , (2.2) 



(2.3) 



(2.4) 



(2.5) 



Note that, as discussed in [11], T)G{n, Nf) can be written as a determinant of a Hankel matrix. 

We regard the Hilbert series and hence T)G{n^Nf) as a partition function of the sys- 
tem consisting of gauge invariant quantities parametrising the moduli space. Let us rewrite 
VG{n,Nf) as 



and interpret H as the Hamiltonian or the energy of the system given by 

n n 

H [T , x) = —2nr'^ log \xa — Xb\ — logwcixc) + n~'^Nf log \T 



(2.6) 



l<a<b<n 



c=l 



c=l 
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n-2 



^ log\xa-Xb\ -^logwcixc) + Nf^log\T ■ 

l<a^b<n c=l c=l 



(2.7) 



The energy of this system can be identified as the potential energy of the two-dimensional 
Coulomb gas, subject to certain external interactions. Consider identical charged particles, 
each with charge +1, located at the position xi, . . . , within an interval (—1, 1) on the x-axis 
in a two-dimensional space. The physical interpretation of each term in the square bracket is 
as follows: 

• The first term is the Coulomb interaction in two spatial dimensions between the particles 
within the interval (—1, 1). 

• The third term denotes the total Coulomb interaction between the charge particle at 
X = Xc, with Xc € (— 1, 1)) and an external particle with the charge —Nf located at the 
position X = T > 1. 

• The second term indicates the total interaction between the charged interval with a 
certain external field determined by the Jacobi weight wg{x). Using the explicit form 
(2.4), this external field is sourced by the extra particles located at a; = ±1 with the 
charges a and /3. 

2.1 The continuum limit 

Let us consider the continuum limit re — )• cxd. We approximate H by using the connection to 
the continuum Coulomb gas (see the books [46, 47]) as follows: 



H ^ % := — - I a{x)\og{T — x)(lx — nT^ / a ix) log WG(x)dx 

n Ja Ja 
/•b i-b 

— dx dx' a{x)log\x — x'\a{x') , (2.1 

J a J a 



where (t(x) is the density of the charged particles supported on the interval a < x < b, 
satisfying the normalisation condition 

b 

a{x)dx = 1 , (2.9) 

and the integration limits a and b are such that 

-l<a<6<l. (2.10) 

Because the (discrete) Coulomb gas contains only positive charges, a{x) is constrained to be 
non-negative throughout the interval [a,b]. 
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In the end we will be interested in the Veneziano limit, i.e., we will take n, Nj — )• oo with 
their ratio kept fixed and of order 1. In view of this, and in order to take into account the 
factor of two appearing for Sp{n) in the integral Do of (2.1), it is convenient to define 



n 



n 
n 



f0TG = Bn,Dn , 

for G = C„ . 



(2.11) 



We can thus write (2.8) for each gauge group as 

% = j dx a{x) v(x) — / \og\x — y\a{y)(iy , (2-12) 

J a \_ J a 

where v(x) is the external potential given by 

v(x) = TZ log(T — x) — log wg{x) 

= 7^ log(r-x)-n~^[alog(l-x)+/31og(l + x)] . (2.13) 

2.1.1 The equilibrium condition 

In order to determine the configuration with the minimum energy (i.e., the equilibrium con- 
dition), we combine the following variation 







/ dx 


\i{x) -2 j 


la 


J a 



5a{x) 



with the variation of the constraint (2.9) multiplied by a Lagrange multiplier —A, 



-A I dx 5a{x) = , A 

J a 



constant 



and equate to zero: 







I dx y{x) — 2 j \og\x — y\a{y)dy — A 

J a \_ J a 



5a{x) 



Therefore, we arrive at the equilibrium condition 



v(2;) — 2 / log \x — y\a{y)dy = constant = A . 

J a 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



This relation indicates that the electrostatic potential, given as the sum of the internal 
Coulomb potential due to the interacting charges within the interval (a, b) and the exter- 
nal potential, is constant at the equilibrium for any point x G (a, b). 

We shall henceforth assume that the system is always at the equilibrium. To avoid 
cumbersome notation, let us slightly make abuse of notation by denoting by Ti the energy of 
the system at the equilibrium. From (2.12), Ti satisfies 



where 



2n = A + I, 



dx v{x)a{x) 



(2.18) 
(2.19) 
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2.1.2 The Hilbert series at large n 

The partition function Vcin^Nf) of (2.6), and consequently the Hilbert series through (2.1), 
can be approximated in terms of the equilibrium Hamiltonian (see, e.g., [45, 48, 49]): 

pG(ri,A^/) ~ e""'^ . (2.20) 

We can thus interpret n^T-i as the free energy of the system, 

F ■.= n^n = ]^n'^{A + I) , (2.21) 

which satisfies 

F \ogVG{n,Nf) . (2.22) 

It is important to control the size of the corrections to the saddle-point result in (2.22). 
Earlier analytic and numerical analyses of similar systems suggest that this result is accurate 
up to, and including, leading and next-to-leading^ order terms in the 't Hooft (large n) [45] 
and Veneziano [12] limits. It is possible to check this numerically by comparing our results 
for the free energy in Section 3 to the exact results in terms of Hankel determinants given 
in [11]. 

3 Exact results in the continuum limit 

Having obtained the equilibrium condition (2.17), we proceed further in Section 3.1 by solving 
this condition for the charge density a in an exact way. Subsequently, in Section 3.3, we use 
this exact expression and (2.21) to obtain the expression for the free energy F. 

As we shall see below, the exact expressions for a and F are functions of the upper and 
lower limits a and h. The latter are determined by degree six equations. In Section 4, we find 
perturbative solutions to such equations in the Veneziano limit. 

3.1 The charge density of the Coulomb gas 

The external force acting on the charged interval can be calculated by differentiating (2.17) 
with respect to x. This yields 

y'{x)-2V r^^dy = , (3.1) 
Ja x-y 

where V denotes the principal value. We aim to compute the charge density cr(x) as a solution 
to this integral equation, with the normalisation according to (2.9).^ This can be achieved 
using the knowledge from the theory of singular integral equations.^ 

^In order for the relation (2.22) to hold at next-to-leading order, it is essential that we are working at a 
critical temperature of the Coulomb gas. At noncritical temperatures, there are finite (but easily calculable) 
corrections even at next-to-leading order (see, e.g., Appendix A in [50]). 

^An alternative method of solving the equilibrium configuration is discussed in Appendix B. 

*See, e.g., Eqs. (25) and (26) in §26 of [51] and Eq. (89.16) of [52] for further details. 
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We follow closely the analysis of [45]. First taking a,/3 > 0, we expect that for the 
unique physical, positive definite solution, the density is supported in the interval (a, b) with 
— l<a<6<l and satisfies the boundary conditions 



or X 



cr{x) — )• as x — )• a 
The solution subject to these boundary conditions is^ 

Vib-x)ix-a) /■S'(x) - v'(y) 



o-(x) 



27r2 



dy 



x-y ^(b-y)iy-a) 



(3.2) 



(3.3) 



where the limits a and b satisfy 



x v'{x) 



dx 



dx 



la y^{b-x){x-a) 2n 
Let us state some comments on the Ansatz (3.3): 



a v^(6 - x){x - a) 27r 



(3.4) 



1. The density (3.3) is supported on the interval (a, b) such that — 1 < a < 6 < 1. They are 
determined by the degree six equations described subsequently in Section 3.2. There, 
we explicitly show that a and b are functions of n. 

2. As discussed in [45], the solution (3.3) can be analytically continued to negative values 
of a and (3, including the values a = —1/2 or /? = —1/2 required for the Bn and Dn 
gauge groups.^ 

Integrating (3.3) to find <t{x) explicitly 
Using (2.13) we obtain 



v'(x) - v'{y) 
x-y 



n 1 

{T-x){T-y)^n 



a 



+ 



{l-x)il-y) (l + x)(l + y) 



(3.5) 



where a and (3 depend on the group G according to (2.5). Then, 



a{x) 



y/{b - x){x - a) 



+ 



n 



27r2 
b 



n 



dy 



a {y-T)^{b-y){y-a) 



a 



+ 



l-x){l-y) (1 + x)(l + y) i ^{b-y){y-a 



dy 



(3.6) 



^The analysis in [45] corresponds to our computation with TZ — Q. The contribution from a non-zero TZ in 
our analysis comes in as the first term of the external potential (2.13). Observe that this term is non-singular 
for all X G [—1, 1], since T > 1. Thus, we indeed expect that the presentation of the density (3.3) of [45] works 
here (as long as q, /3 > 0). 

®The analytic continuation works up to high-order corrections in the Veneziano limit. It would also be 
possible to start with negative a and P, and analytically continue the result to the positive values. In this case 
the Ansatze (5.7) and (5.13) in Section 5 would be relevant. 
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Using the integral (A.l), we find that the density a{x) is given by 



2TTa{x) 
{b - x){x - a) 



1 

n 



a 



+ 



7^ 



(r-x)V(T-a)(r-6) ' 

Let us make a few comments on (3.7): 



X G (a, 6) . 



(3.7) 



1. The terms in the first line of (3.7) are equal to the density computed in [45], which 
correspond to our analysis with TZ = 0. The contribution from a non-zero TZ comes in 
as the second line of (3.7). Observe that the latter does not depend on a and (3. 

2. Recall that we take first a,/3 > 0. Since — l<a<a;<6<l and T > 1, the second 
term of (3.7) is strictly negative, whereas the first term (coming from the Jacobi weight) 
is strictly positive. In order to obtain a positive charge density in the large n limit, it 
is crucial that 1) the first term cannot be neglected and 2) the limits a and b must also 
be certain functions of n in such a way that the terms in the first line can balance that 
in the second line. We compute a and b perturbatively in the subsequent section; see 
Eqs. (4.2) and (4.14). 

3. As pointed out in [45], the partition function Vcin^Nf) computed using (3.7) can be 
continued analytically to negative a and /3, due to the real analyticity of the expres- 
sion. For convenience, we shall henceforth assume that a, /3 > and, at the end, we 
analytically continue the result to all cases listed in (2.5). 

3.2 The integration limits in the Coulomb gas 

Let us determine the constraints on the integration limits a and b. The first equation of (3.4) 
reads 



27r = 7^ 



X dx 



a 



X dx 



a {x-T)y'{b-x){x-a) nJa {x - l)y^{b - x){x - a) 
13 xdx 
nJa {x + l)^{b- x){x - a) ' 



(3.8) 



Performing partial fractions and using identity (A.l), we obtain 



a + B _ 

2 + — --n + 



TIT 



n 



^J{T-a){T-b) 
The second equation of (3.4) reads 
„ . dx 

o = n 



1 

n 



a 



+ 



/3 



V(l-a)(l-6) ^(l + a)(l + 6) 



(3.9) 



a 



a {x-T)y^{b-x){x-a) nj^ (x - 1) v^(6 - x)(a; - a) 



dx 
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dx 



la {X + l)^y{b- X){x - a) 

Using identity (A.l), we obtain 



(3.10) 



n 



ViT-a){T-b) 



1 

n 



a 



V(l-a)(l-6) V(l + a)(l + 6) 



(3.11) 



Substituting (3.11) into (3.9), we obtain 



a + /3 _ 
2 + -TZ 



n 



1 

n 



a{l-T) 



+ 



/3(i + r) 



V(l-a)(l-6) y/(l + a)(l + 6) 



(3.12) 



In subsequent computations, we make use of (3.11) and (3.12) as two constraints on the 
integration limits. Note that these constraints are indeed the generahsation of the result in 
[45] to the case of non-zero TZ. 



Change of variables To proceed further, let 

X:= Y 



^(1 -")(!- i") ' 

and express ah and a + 6 in terms of X and Y , 



V(l + a)(l + 6) 



ah 



- — — 

2 V X2 + y2 



Using (3.11) and (3.14), we obtain 

^2 1 / 1 1 



T + 



1 



a + h 



1 1 

+ 



1 / 1 



2 V y2 x2 



1 



n27^2 



(3.13) 



(3.14) 



(3.15) 



2 x'^j 2 \x'^ y2 j ^f^x - /3y)2 ■ 

From (3.12) and (3.13), we also have a /inear relation between X and Y: 

-a{T - l)X + P{T + 1)Y = -n{n - 2) + (a + /3) . (3.16) 
We shall make use of (3.15) and (3.16) in subsequent computations. 
3.2.1 The degree six equation 

Solving for X from (3.16) and substituting it into (3.15), we find that, upon disregarding the 
solution T = 1, y satisfies the following degree six equation: 



(3.17) 



fc=0 



where the coefficients are given as follows: 

C6 = 8/3^f 3 , 
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C5 



33r7-i2 



T^a + P -Tin] (2 + T 



C4 = Alf]- nTZ) (2a + /3- nTZ) + 2T\ W + (a + /3)(5a + 3/3) 



+ n7^ (^4(n - 2(a + /?)) + 5n7^) | + (2n + a + /3) (-2n + a + /3 - 2n7^) 



C3 = -8/3r a + /3-n7^ 



C2 



4n^ + 2a/3 + 2n7^ ( 2n - a - /3 + n7^ 



+ ( -2n^ + a{a + /3) - n(2n + a + /3)7^ ) T 
2 r 



a + /3 — nTZ 
+ 2n(2n + a + 



16n2 - 4/3^ + 4n27^ (^4 + 7^) - 2Tf^An'^ + /3(-2a + 3/3) 
/3)^} + (a-/3)(a + /3)r2^ 



Cl = 2/3 a + /3-n7l 2 + r 



Co 



a + /3 — 



with 



n = n-2 



T + l . 



(3.18) 



(3.19) 



One can also similarly compute the degree six equation for X. 

We conjecture that these degree six equations arise from the existence of the Painleve VI 
equation that the Hilbert series satisfies [11]. 

3.3 The free energy 

Recall from (2.21) that the free energy receives two contributions, from A and /. 
Computing A from (2.17) 

We compute A as follows. First, let us multiply both sides of (2.17) by the factor ^[(6 — 
x){x — a)]^^/^ and integrate over x: 

=dx 



A 



dx 



vr J a y^(6 - x){x - a) 



v(x) 



a TT^ib - x){x - a) 
^ dx 



TTy^ {b — x){x — a) 



I log|x - y|cr(y)dy , (3.20) 

J a 



and so 



A 



f- 

J a vr 



v(x) 



dx-2 I dy a{y) 



log\x - y\ 



a VTy^ {b — x){x — a) 



dx 



(3.21) 



Note the following important observation (see, e.g., (6.19) of [53]): 



d 



log \x — y\ 



dy Ja y^{b-x){x - a) 



dx = V 



dx 



a (y - x)^/{b - x){x - a) 



(3.22) 
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if y € [a, b]. In other words, the integral with respect to x in the second term on the RHS of 
(3.21) is independent of y. Hence, we can replace y in the logarithm by a constant, say a: 



log|x-y| 



zdx 



a •y/(6 - x){x - a) J a y/{b - x){x - a) 



log |a; — a| 



zdx 



(3.23) 



Therefore, 
A-- 



" v{x) 
vry^ {b — x){x — a) 

^ v(x) 

^T^y {b — x){x — a) 
^ v(x) 



dx - 2 



dx - 2 



dy a{y) 



log{x — a) 



a 7Ty^ (6 - x){x - o) 



a ^T^y{b- x){x - a) ^ 
dx 



TT 



yj{b- x){x - a) 



dx — 4 log ( — a 



(3.24) 



Observe that this expression of A is independent of the charge density a{x). 
Using (2.13) and (A. 8), the integral in the first term is 



J a n 



v(x) 



y/{b- x){x - a) 



dx = 27^1og 



2 

n 



a log 



+ /31og 



yi + g + yi + b 

2 



(3.25) 



Thus, we obtain 



A = 2n log 



4 log ( -Vb^ a 



2 
n 



a log 



+ /31og 



(3.26) 



Computing the integral / in (2.19) 



dx y{x)a{x) 

= n ( dx log(r - x)(j{x) -- [ dxa(x) [a log(l - x) + ^ log(l + x) 

J a n J a 

Using (3.7), we have 

I3n 



(3.27) 



2TTn 



an 



v/(l-«)(l-&) ^{l + a){l + h) n^{\-a){l-h) n^{l + a){l + h) 

a/3 I a/? 



iV(l-a)(l-&) ny/{l + a){l + by 
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n 



I,'R--Il^^-It 
n n 



(3.28) 



2^^{T-a){T-b) 
where the integrals if , if^ ,Hilf are defined as 

T± _ cb ^^ '^osiT-x)y/ib-x){x-a) ± _ ^^ iogi'i--^)Vib-x){x-a) 

-'l - Ja {l±x) ' ^2 - Ja '^^ l±x 

T+± — rb (b-x)(x-a) _ ^^ iogiT-x)y/{b-x)(x-a) (3 29) 

^2 - Ja l±x ' - Ja (T-x) ^ ' ' 

T± _ j-b ^^ log(l±x)-^ {h-x){x-a) _ 

^4 - Ja '^^ T-x ' 

they are computed explicitly in Appendix A. 2. 

4 Series expansions in the Veneziano limit 

In this section the main goal is to compute, from the preceding exact results, the charge 
density, the free energy, and hence the Hilbert series perturbatively in the Veneziano limit, 
i.e., as a power series in 1/n while keeping the ratio TZ defined in (2.11) fixed and finite. 

As it turns out (see also the analysis of Section 5), the result is essentially different 
depending on whether the system is in one of the following phases: 

1. The gapless phase: The name "gapless" denotes the phase in which the integration 
limits a and h are such that a — >■ —1, 6 — )• 1 as n — )• oo. In this phase, the charge density 
is a continuous function supported on the interval (—1, 1) in the large n limit. In other 
words, there is no gap created on the charged interval even in the presence of external 
interactions. 

2. The gapped phase: On the other hand, in the "gapped phase", the lower limit a 
tends to oq, for some oq strictly greater than — 1, as n — )• oo. In the latter phase, the 
gap is created in the interval (— l,ao) and the charge density is supported only on the 
interval (oq, 1) in the large n limit. The presence of such a gap is due to the attractive 
interaction due to the external particle at x = T. 

Subsequently, we show that the relevant parameters involved in these two phases are IZ 
and T. Explicitly, we have the following: 

• The gapless phase: < 7^ < 2, or 7^ > 2 and T > Tc (i.e., < t < Q , 

• The gapped phase: TZ > 2 and 1 < T < Tc {i.e., tc < t < I) , 
where the relation between T and t is given in (2.3) and 

t-^ ^_^(^^2) + 2 

^'^"Te-i' '~ 2(7^-l) • ^ ' 

Notice that taking TZ = 2Nf/n for the gauge groups U{n) and SU{n), the critical value tc 
agrees with that of [12]. 

We emphasise again that in the following computation a, /3 > is assumed; the end result 
is then analytically continued also to the cases a,/3 = —1/2. 
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4.1 The gapless phase 

The first task is to find the endpomts a and b perturbatively at large n. In this phase, we 
require that a ^ —1 and 6 — )• 1 as n — )• oo. We therefore propose the following Ansatze: 

« = -l + 55 + ^5 + o(i_) , 6=i_»|_1+o(i_) , (4,2) 

where 02,03,62,63 are fixed and finite such that 02,62 > 0, so that the interval [a,b] is 
contained in [— 1, 1] as n — )• 00. 

Inserting these Ansatze in (3.15) and (3.16) and recalling the definitions (3.13), we find 
the perturbative solution 



(4.3) 



_ 2(T+l)2/32 _ 2(T-l)2a2 

~ [i2-R)iT+l)-RVT^]^ ' ^ ~ [(2-R){T-l)+RVT^f ' 

^ 4(T+l)3/3^(«+/3) ^ MT-Ifa\a+I3) 

^ [(2-R)(T+l)-Rv/T2^]^ ' ^ [{2-R){T-l)+RVT^f 

Observe that in this solution, 02 and 62 are strictly positive for T > 1, as required. 

4.1.1 The charge density 

Substituting (4.2) into (3.7), we obtain 

for all X £ (—1,1). We analytically continue this result to all cases listed in (2.5). 

Observe that the leading contribution does not depend on a and /3. It is non- negative 
for all X G [—1, 1] if and only if 

TZ(TZ — 2) + 2 

< 7^ < 2 or 7^ > 2 and T > Tc := — '-y- ■ (4-5) 

When the subleading contribution is included, we consider each follows. 

• For G = Cn, a = P = 1/2 and so the subleading term is strictly positive."^ 

• For G = Bn, a = —f3 = 1/2 and so the subleading term as written in (4.4) vanishes. 
Formally, the subdominant contribution suppressed by 1/n^ can be computed by ex- 
panding a and 6 further to include terms at order O {n'^) ■ However, the result for the 
latter is too lengthy to be reported here, and it is not likely that the result at such 
an order will be useful, because the relation (2.22) is expected to receive nontrivial 
correction at high orders. 



'^This statement is true for all x within the open interval (—1, 1). The expansion is involved at the endpoints, 
which is reflected by the fact that (4.4) does not satisfy the normalization condition at next-to-leading order. 
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For G = Dn, a = P = —1/2 and so the subleading term, suppressed by 1/n, is strictly 
negative. 

One can ask when the perturbative density (4.4) breaks down; in other words, when 
cj(x) contains a zero for some x E [—1,1]. The zero of cr{x) occurs at 



2-n 



+ 



n 



7^^/^2 - 1 



(4.6) 



For < 7^ < 2, it is immediate that Xc > 1, which hes outside the support [—1, 1] of the 
density; thus, a{x) is strictly positive for all x G [a,b]. On the other hand, for 7^ > 2 
and T > Tc, Xc G [—1)1] and the charge density thus becomes negative if and only if 

27^(7^-l) + l , . ^ . (7^-2)(l + T) + 7^^/^^ 

ir < < :^ and 1 < n < , — ^ —^^ — . (4.7 

27^-l - 2[2(7i- i)(r + i) -7^2] ^ ^ 

Since n is bounded from above by a finite number, higher orders in the perturbation 
theory become significant. 

These arguments support the validity of analytic continuation we performed earlier. 
4.1.2 The free energy 

Let us now compute the free energy F given by (2.21). There are two contributions, from the 
constant A and the integral /. 



The contribution from A The constant A of (3.26) is 

^ = [2 log 2 - 7^ log(2t)] + - (a + /3) log 2 + 4^ A^"^^ + O (n~^) 



n 



where t and T are related by (2.3) and 



A 



(-2) 

A 



+ 



(l + t)/32 



(4.8) 



(4.9) 



(7^-l)^ + l (7^-l)^-l 

It should be observed from (4.8) that the leading term does not depend on a and /3 and is 
hence universal for all groups listed in (2.5). 

The contribution from / The integral / is defined in (3.27). A direct computation yields 
I = n[n log (1 - t^) - log(2t)] + ^ [(a + /?) log 2 

- 27^{a log(l - t) + /3 log(l + t)}j +C'(n"2) . (4.10) 
The free energy By inserting these expressions in (2.21) the free energy becomes 



gapless 



n 



log2-7^log(2^) + -7^2log(l-^2) 



+ n 



[a + 13) log 2 



n{a log(l -t) + l3 log(l + t)}]+0 (nO) 



(4.11) 



Again, we emphasise that the leading contribution does not depend on a and (3 and is hence 
universal for all groups listed in (2.5). 
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4.1.3 The Hilbert series 

We can now readily compute the partition functions and hence the Hilbert series from (2.22). 

1. 



log Vg ~ --Fgaple 



-n 



log2-7^log(2^) + ^7^2log(l-^2) 



n 



{a + p) log 2 

-7^{alog(l-^) + /31og(l + ^)}J +0(n°) . (4.12) 
Using (2.5), (2.11) and (2.1), we find that the Hilbert series are given by 

' '\Nf{Nf + 1) log(l - t2) , G = Bn, Dn 



log gN!,G{t) 



-{2Nf-l)Nf\og{l-t') , 



G — Cn 



(4.13) 



Observe that these asymptotic results turn out to coincide with the exact results given by 
[7, 11] for the case < 7^ < 2. In addition, the negative coefficient of log(l — t^) indeed 
coincides with the dimension of the moduli space for any 7^ > 0. 

4.2 The gapped phase 

In the gapped phase, we take Ansatze for a and b to be 



a = ao + 



ai 



O 



n 



1 



(4.14) 



where 00,01,62,63 are fixed such that oq > — 1 and 62 > 0, so that the interval [a, 6] is 
contained in [—1,1] as n — )• 00. 

These coefficients are determined by substituting these Ansatze to (3.15) and (3.16). In 
order to calculate the charge density and free energy at leading order, only the coefficients 



ao 



IZ"^ - 4:(n - l)T 



(T- l)a2 



(4.15) 



(7^-2)2 ' " 4(7^-l) 

are needed. Notice that oq is the same as in (B.19). For the subleading order, it is necessary 
to take the following coefficients into account: 

27^2(^ - l)[S{a + (3)-{n- 2)/3] 



ai = -- 



63 



(7^-2)35 

a^(r - l)[(7^ - 2)a + (7^ - 5 - 2) 
8(^^0)2 



where we define 



S = y^23^(^3I)(r:^r) = ^/2{n-l){Tc-T) . 



(4.16) 
(4.17) 

(4.18) 



Note that the quantity in the square root in (4.18) is non-negative, ao > — 1 and 62 > 
if and only if 

(1< T < Tc and 7^ > 2) or (1< T < and 1< 7^ < 2). (4.19) 
In the following, we use the positivity of the charge density to rule out the second condition. 
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4.2.1 The charge density 

Using the above coefficients, we can compute the density from (3.7) which yields 

where ag is given in (4.15) and the subleading term ^\x), suppressed by 1/n, is given by 

= (^-2)(l + x)(a + /3)-2g/3 
27^(7^ - 2)(1 + x) V(x - oo)(l - x) ' 

Observe that the leading order of o"(x) is positive definite if the first condition in (4.19) 
holds. Hence, the condition for the gapped phase is 

KT <Tc and 7^ > 2 . (4.22) 

We shall henceforth assume this condition throughout this section. 

Note that the charge density as given in (4.20) can be analytically continued to all cases 
listed in (2.5). 

4.2.2 The free energy 

Let us now compute the free energy. There are two contributions from the constant A and 
the integral /. 

The potential term A of (3.26) is 

A = 2{n- 1) log(7^ - 1) - 2(7^ - 2) log(7^ - 2) + (7^ - 2) iog(r - 1) 

+ Ia^-i) + O f i") , (4.23) 

and the subleading term, suppressed by 1/n, is 

A^"^^ = /? log 2 + 2(a + /3) log(7^ - 2) 

-2/31og(7^ + 5-2)-alog{(7^-l)(^-l)} . (4.24) 

The integral I defined in (3.27) becomes 



-n 



/(7^-l)2^ ^ /(7^_l)^T^ 



+ Ia;-!) + O (n-2) , (4.25) 
n 



and the subleading term A^ ^■^ is 

A^~^) = -{n - l)/31og 2 - 2(7^ - l)(a + /?) log(7^ - 2) + 27^(a + log(7^ - 1) 
+ 2(7^ - l)/3 log(7^ + S-2)-a log{(7^ - l)(r - 1)} 

-7^/3log{l + 7^(7^ + 5-^-l) + ^} . (4.26) 
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Substituting the expressions for A and / in (2.21), the free energy is 
1 



gapped 2 ^ 



log{7^(7^ - 2)} + 27^log(T - 1) - 2 log{(7^ - i)(r - 1)} 
- 4(7^ - 1) log(7^ - 2) - 2(7^ - 2)7^log(7^ - 1) 

-n (7^-2)/31og2 + 2(a + /3)(7^-2)log(7^-2) - 27^(a + /3) log(7^ - 1) 
2(7^ - 2)/31og(7l + S - 2) + 2alog{(7^ - l)(r - 1)} 



+ n/3 iog{i + n(n + s-t-i) + t} 



(4.27) 



The leading order term is again independent of a and /3 and thus universal. Interestingly, 
one can check that this term (as well as the one in gapless phase, (4.11)) also agrees with the 
one for the U{n) and SU{n) gauge groups (given in Eq. (55) of [12]), after the leading order 
factors from the relations (2.1) are included. We also emphasise that the free energy (4.27) 
can be analytically continued to all cases listed in (2.5). 

4.2.3 The Hilbert series 

We use (2.22) to compute the partition function and the Hilbert series follows easily from 
(2.1). For illustration, let us present the result in the limit t — )■ 1, where the expressions 
simplify; we set 



t = l- r] 



with r/ — 7- 



(4.28) 



We thus obtain 



log gNf,c„ = nlog(2) + [2^^ - + 2(^ " log(^ - 2) 



+ 



1 



(7^ - l)n + (7^ - l)^n2 log(7^ - 1) 



- -n^n^ log{{n - 2)7^} - ^2(7^ - l)n'^ - n] log(?7) + O {rj; n°) , (4.29) 



log gN 



-n log(2) + ^ (7^ - 2)n + 2(7^ - l)n^ log(7l - 2) 



+ 



1 



{lZ-l)n + {lZ-lfn^ log(7^-l 



loggNf,D^ 



n^r?\og{{n-2)n} - n{n-l){2n+l) log(r/) + O (r?; n°) 
1 



(4.30) 



+ 



nlog(2) + 
rl 



-(7^-2)n + 2(7^- l)n^ log(7^ - 2) 



1 



{IZ - l)n + (7^ - lyn^ \og{n - 1) 
n^n^ log{(7^ - 2)7^} - \2{n - l)n^ + n] log(?7) + O (r/; n°) 



(4.31) 
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Note that the minus of the coefficient of logry = log(l — t) in each case indeed coincides with 
the dimension of the moduh space as given in [7, 11]. 

Moreover, we can perform further non-trivial checks of these expressions by considering 
the following examples: 

• G = Cn with n = 2 + 2n-^ ov Nj = n + l. Eqs. (4.29) and (4.27) of [11] are in 
agreement with the limits — )• 0, n — )• oo: 

log5^,,c„ = log(l - t^^n - Nf{2Nf - 1) log(l - t") 

n(3 + 2n) log(2r?) . (4.32) 

• G = Bn with 7^ = 2 + or Nf = 2n + 1. Eqs. (4.30) and (2.36) of [11] are in 
agreement with the limits — )• 0, n — )• oo: 

log <7iV,,B„ = log(l + t^O - \Nf{Nf + 1) log(l - t") 

~ -n(3 + 2n) log(27?) . (4.33) 

m G = Dn with IZ = 2 + n-^ or Nf = 2n + 1. Eqs. (4.31) and (3.31) of [11] are in 
agreement with the limits 77 —t- 0, n — )• 00: 

log ff7V,,B„ = log[l + (1 + 2n)t2« - (1 + 2n)t2"+2 - t4n+2] _ ]^Nf{Nf + 1) log(l - 

~ -n(3 + 2n) log(2r/) . (4.34) 



5 Leading order approach in the Veneziano hmit 

In this section, we present an alternative approach to derive most of the results to the Hilbert 
series presented above. The purpose is to bypass the exact results presented in Section 3 and 
inferring the results more simply by using standard perturbation theory; see [40, 41, 50, 54]. 
Solving the equilibrium condition at leading order (LO) in the Veneziano limit is enough to 
obtain the free energy and hence the Hilbert series up to next-to- leading order (NLO), as we 
will demonstrate below. 

To be precise, we write the potential as 

v(x) =vo(x) + -A^'^(x) , (5.1) 
n 

where 

vo(x) =7eiog(r-x) (5.2) 
is the leading contribution in the Veneziano limit, and the contribution from the Jacobi weight 
a1 (x) = — log wg = —a log(l — x) — (3 log(l + x) (5.3) 
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appears at NLO. Then, generalizing the notation of the previous sections, we can expand the 
charge density 

a{x) = ao{x) + -A(-i)(x) + ^A^'^Hx) + ■■■ . (5.4) 

Below we will use analogous notation for series coefficients for other quantities like A and /. 
The terms in this expansion can then be solved by applying standard perturbation theory to 
the equilibrium condition (2.17). However, as we already mentioned above, for our purposes 
it is enough to restrict to the solution at LO.^ As is usual in perturbation theory, the NLO 
correction to free energy can be obtained in terms of the LO solution and the NLO "source" 
of perturbation, i.e, the Jacobi weight of (5.3). 

5.1 Leading order solution 

In the LO computation, only the LO piece (5.2) contributes to the external potential. The 
equilibrium condition (2.17) thus reads at LO 

vo(x) - 2 / dy log \x - y\ao{y) = const. = , (5.5) 

J a 

where the potential is given in (5.2), and the support of cro{x) is [a,b]. 

Before going to the explicit solutions to (5.5), let us discuss the physical interpretation 
of the system. Recall that cr{x) contains only carriers with charges having a single sign, and 
therefore it cannot become negative. The external charge at x = T > 1 has the opposite 
charge to the charge density o'(x) on the interval —1 < x < 1, and therefore attracts it. 
Therefore, we expect that we will always have 6 = 1. We further expect to have two kinds of 
solutions to (5.5) depending on the values of T and TZ: 

1. For a sufficiently large T or sufficiently small TZ, the attraction is not sufficiently large 
to move the end point x = — 1 of the interval and hence the lower limit a = — 1 in 
this case. Note that in this case, the charge density (7o{x) blows up at both end points 
X = ±1. 

2. As 7^ is increased and T is decreased to sufficiently large and small values, respectively, 
the attraction drags the end point x = — 1 of the interval to x = ag > —1 and therefore 
a gap appears for x G [— l,ao). When the gap has appeared, the charge density must 
vanish continuously at x = oq (but it still blows up at x = 1), since otherwise there 
would be a net force acting on the charges near the edge of the gap, so that moving the 
charges into the gap would lower the energy of the system. 

The first situation is known as the gapless phase and the second is known as the gapped phase, 
which were already discussed in Section 4. We will determine oq and the values of TZ and T 
at the transition point below. 

^Finding the perturbative solutions at high orders is involved, in particular when a and /3 are positive, due 
to singular behavior near the endpoints of integration. At our LO analysis such complications are absent. 
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5.1.1 Charge density in the gapless phase 

According to the above discussion, we expect that for the unique physical, positive definite 
solution to (5.5) in the gapless phase, the density is supported in the interval [—1,1] with 

f7o(x) — )• +00 as X ^ ±1"^ . (5-6) 

The solution of (3.1) subject to the boundary conditions (5.6) is^ 

Mx) = ^ f ^ yr37dy+ -JL= . (5.7) 

27r2vi - J-iU - X VI - x^ 

Here the constant ci can be determined from the normalisation condition (2.9). Note that 
the integral over x from —1 to 1 of the first term in (5.7) vanishes due to the integral identity 
(A. 2). We therefore obtain 

ci = l/vr . (5.8) 
Substituting (5.2) into (5.7) and using the identity (A. 3), we find 



^o(a;) = „ ^ , V \ 7 77 ^dy + 



27r2^/r^ y_i (y-x)(y-r) ^ ^VT 

(2 - 7^)(^ - x) + 7^VT2 - 1 



X 



2 



27r(r - x)\/r 



X 



2 



(5.9) 



Conditions for gap formation Let us then discuss when the solution (5.9) is valid. If 
7^ < 2, it is easy to see that (y{x) is nonnegative for all acceptable values of t or T . Therefore 
it is the desired solution. For TZ > 2, the numerator of a takes its smallest value at x = —1. 
Requiring this to be nonnegative leads to 

(2 - 7^)(^ + 1) + 7^\/T2 - 1 > o or ^^>n. (5.10) 

This translates to 

2(7^- 1) - - n-i ^ ' 

5.1.2 Charge density in the gapped phase 

We expect that the physically reasonable solution in the gapped phase satisfies 

ao{x) — )• +00 as x — )• 1~, and (7o{x) — )• as x — )• a^}' . (5-12) 
The unique solution of (3.1) subject to the boundary conditions (5.12) takes the form^*^ 



1 Ix-ao /"^ v'o(y) / 1 - y 



ao 



y-x y y -ao 



^See, e.g., Eq. (26) on Page 131 of [51]. 
^°See, e.g., Eq. (25) on Page 131 of [51]. 
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In order to compute oq, we use the normalisation condition 



1 = / cJo(x)d2; = I yo{x)J- — -dy = ^ [ ^f^^J- — -dy 



Thus, we obtain 

ao = (^_2)2 • (5-15) 
Substituting (5.2) into (5.13), we obtain the leading term of the density to be 

= T3^^ Ao (y-T)(,-x)V^'^ 



(7^-2)v^^ 



27r{T - x)y/T^ ' 

where the last equality follows from partial fractions and identities (A. 4) and (A. 5). 



(5.16) 



Conditions for gap formation In the gapped phase, (5.16) is positive when TZ > 2. The 
condition — 1 < ao < 1 for the expression (5.15) further leads to 

KT <Tc or tc<t<l. (5.17) 

Notice that these conditions are complementary to those of the gapless phase above. We 
have thus identified the LO solution for all acceptable parameter values. 

5.1.3 The free energy 

In order to compute the LO terms for the free energy in (2.21) we need to compute A of (2.17) 
and / of (2.19) at LO in the Veneziano limit. 

Gapless phase 

The LO result for A can be calculated from its definition (5.5). The calculation can be 
simplified by using the identity (3.24), which, by its derivation, holds for any potential v(x) 
and charge density a having the support [a, b]. Applying the identity to the LO equation (5.5) 
with a = — 1 and 6 = 1, we obtain 



-TZ log(2r - 2\/r2 - 1) + 2 log 2 = -7^ log(2t) + 2 log 2 . (5.18) 



^'^In the critical case TZ > 2 and t — tc, we found that both sohitions (5.9) and (5.16) are valid, and in this 
case the two expressions indeed match. 
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The LO integral Iq is computed directly from (2.19), giving 



vo{x)aoix) dx 



(TZ - 2)n 



log l2T - 2VW^l) - log (^L= + 



= -7^ log (2t)+7^2 log (1-^2) , 

where we used the integral identities (A. 8) and (A. 14). 
Gapped phase 

Setting a = aQ and 6 = 1 in (3.24), we obtain 



(5.19) 



An 



n 



ao 



log(T - x) 
x){x - ao) 



drc — 4 log 



= 2(7^ - 1) iog(7e - 1) - 2(7^ - 2) log(7^ - 2) + (7^ - 2) iog(r - 1) 

= 7Jlog(R-l) + (R-2)log[(^^^iiil=f' . 
The integral Iq is computed as in the gapless phase by using the definition (2.19): 



(5.20) 



/ vo(x)cro(x) 

J ao 



dx 



7^(7^-2)v^-ao 



log 



1 



1 



2^T - ao 2VTTT 



7^(7^ - 2) log 



ao 



(l-t)'^ lZ-2 

7^(27^ - 1) log 2 + 7^ log ^ ' + 7^(7^ - 2) log — — 27^(7^ - 1) log 



4t 



5.2 Linear perturbation to the free energy 



7^ 



2 

271-2 

n 



(5.21) 



Let us then discuss how the NLO correction to the free energy can be obtained. Expand- 
ing (2.12) up to NLO we find that (see also, e.g., [49]) 



n = ^{Ao + lo) + -a57^) + O 



n 



1 



(5.22) 



where 



A, 



(-1) 

H 



/ ao{x)Ai-^\x)dx 

J a 



a / (To(x) log(l — x) da; — /3 / cro(x) log(l + x) dx . 



(5.23) 



Here a equals —1 in the gapless phase and oq, as given by (5.15), in the gapped phase. Notice 
that this expression is independent of the NLO solution "^^(x). This was expected: o"o(x) 
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extremizes the Hamiltonian at LO, and therefore any variation around it leaves the Hamilto- 
nian unchanged at hnear order. Therefore, only the LO solution to the equilibrium condition 
is required to compute the variation of the Hamiltonian ^\ as is usual in perturbation 
theory. 

The gapless phase The NLO corrections to the free energy are obtained by direct inte- 
gration of (5.23). We need to evaluate the following integrals with fTo(x) given by the gapless 
phase case of (5.9): 



1 



2^ y_i \^ T-x J y/l^ 

^-ho,2-^U^^^) (5.24) 



2 ^ 2 ^Vt + ^/T^ 
= log2-7^1og(l±^) . 

Therefore, from (5.23), the NLO correction to the free energy is given by 

aJ7^^ = a J- + /3J+ 

= a[log2-7^1og(l-^)] + /3[log2-7^1og(l + ^)] . (5.25) 

The gapped phase In the gapped phase we proceed by direct integration as well. The 
relevant integrals are as follows with o'(x) given by the gapped phase case of (5.16): 



I± := — / a{x) log(l lb x) dx 

n-2 1 jx - ao 



2vr A„(r-.)V l-x 

~ t - il . (5.26) 



ao 



T-x ) ^(l-x)(x-ao) 



The relevant integral identities are (A. 16) and (A. 17). With the aid of these identities, we 
obtain 

2:- = ^^o-|iog(r-i) 



l+ao ^ n V l+ao 



where Aq is given in (5.20). Substituting oq from (5.15) and simplifying, we obtain the 
correction to the free energy from (5.23) as 
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- (7^-2)/31og2 + 2(a + /3)(7^-2)log(7^ 



2) - 27^(a + P) log(7^ - 1) 



2(7^ - 2)/3 log(7^ + 5 - 2) + 2a log{(7^ 



i)(r-i)} 



+ 7^/3 iog{i + n{n + s-t-i) + t} 



(5.28) 



where S = y^TZ'^ - 2(71 - + T). 



5.3 Discussion 



The free energy in each phase can now be obtained by using 



F = -{Ao + lo) + nA^^'> + O (n' 







) 



(5.29) 



and the results in (5.18), (5.19), and (5.25) [in (5.20), (5.21), and (5.28)] for the gapless 
[gapped] phase. 

As a consistency check, all results obtained in this Section match with those of Section 4. 
In particular, 

• The charge densities (5.9) and (5.16) equal the leading order terms of (4.4) and (4.20), 
respectively. 

• The results for the free energy, as obtained from (5.29), match with (4.11) and (4.27) 
in the gapless and gapped phases, respectively. 

Thus, we also see that inserting the free energy of (5.29) in the relation (2.22) reproduces the 
key results of this article, given above in (4.13) and in (4.29)-(4.31). 

6 Phase transition 

Having been discussing the two phases of the moduli space in the preceding sections, we now 
examine the transition between these two phases. Let us assume throughout this section that 
7^ > 2 so that the transition is present. Notice that the large n limit is also required for the 
transition to take place: for any finite n the integral T>g in (2.1) is an analytic function of 
its parameters, and hence no phase transitions are expected. At finite but large n, we thus 
expect that the system is well described by the (leading order) results in the Veneziano limit 
except for very close to the transition point, where the phase transition is smoothed out. We 
can use the next-to-leading corrections, which we have calculated above, to estimate the scale 
where the smooth-out takes place. Let us now discuss in detail how this works. 
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6.1 The gap width near the transition point 

We first analyse the width of the gap in the gapped phase near the transition point. Indeed, 
this width can be regarded as the order parameter of our system: if it vanishes (at leading 
order in the Veneziano limit), the system is in the gapless phase, otherwise the system is in 
the gapped phase. 

Consider the system in the gapped phase. The gap spans the interval [— l,a] and hence 
the gap width is L := a — (—1) = a + 1. Using (4.14), (4.15), and (4.16), we obtain 

T^^j-Mn-1)T 27^2(^-l)[g(a + /3)-(7^-2)/3] 
^^' + ^^^2? ^^^2)35 + 0{n ). (6.1) 

Define 

e := Tc — T , T < Tc in the gapped phase (6-2) 
and consider the limit T approaches Tc from below or e — )• 0"*". We obtain 

7^2/3 1 7^2(a + /3) 



m-ll 1 

J J f -4- — 

(7^-2)2'^n 



^/2(7^ - 1)3/2 (7^-2)(7^-l) 

+ O (n-2) , e ^ 0+ . (6.3) 



We encounter here a double scaling problem, since n and e can be varied independently. We 
can identify two regimes: 

1. The case that n^^e^^/'^ <^ e, or equivalently ne^/^ — )• oo. The first term of (6.3) is 
dominant. The order parameter near the transition point is therefore 

2. The case that e ^ n^^e^^^"^, or equivalently ne^/^ — 0. The second term of (6.3) is 
dominant over the first term, suggesting that the perturbation series (in 1/n) diverges. 
We interpret this as a signal of the phase transition having disappeared, so that the gap 
size L behaves smoothly over the transition region. 

The scale at which the phase transition is smoothed out is thus e ~ n~2/3_ 
6.2 Free energy at the phase transition 

Let us consider the difference between the free energies of the two phases: 

• — -^gapped -^gapless (6'5) 



^ We have verified this numerically for the Coulomb gas system of the Sp{n) gauge group defined in Section 3. 
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in the hmit 



e := 



as n 



oo . 



(6.6) 



Using (4.11) and (4.27), we obtain 
8(7^- 1)5 



6F 



37^4(7^- 2)4 

+ 



^3 _ 12(7^-^lm7^-2) + 2] ^4 ^ ^ ^^5| 



7^6(7^-2) 



6 



n 



8^2(^-1)5^ 3/2 2(^-l)3(a+^ / /2 

37^2(7^-2)2 ^ 7^2(7^-2)3 ^ ^^{^ 



n 



+ O(n0) 



Let us now consider 5F in the cases discussed in the preceding subsection: 
1. The case ne'^/2 _). The leading contribution to (6.7) is given by 



6F 8{n - 1) 



3n^{n - 2Y 



(7^-2)2(7^- 1) 

2^^4 



• 



(6.7) 



(6.8) 



where the order parameter Lq is given in (6.4). This case exhibits the third order phase 
transition, in a similar fashion to [18, 19]. 

2. The case ne^^"^ — )• The terms on the second row of (6.7) are leading with respect to 
the ones on the first row, which again suggest that the perturbative results cannot be 
trusted. We expect that there is no phase transition but a smooth cross-over in this 
limit. 



Acknowledgments 

Y. C, N. J., and M. J. thank the Max Planck Institute for Physics for warm hospitalities while 
this work was in progress. N. M. wishes to express his sincere gratitude to Amihay Hanany for 
numerous discussions since 2008 that triggered his interest in this direction. He also thanks 
the hospitality of Universidade de Santiago de Compostela, CERN Winter School 2013, ICTP 
Spring School 2013, and the Fizzi-Fazzi doublet (Fabio Apruzzi and Marco Fazzi) during his 
visit. The work of N. M. is supported by a research grant of the Max Planck Society. N. J. is 
supported by the MICINN and FEDER (grant FPA2011-22594) and the Spanish Consohder- 
Ingenio 2010 Programme CPAN (CSD2007-00042). N.J. is also supported by the Juan de 
la Cierva program. The work of M.J. was in part supported by grants PERG07-GA-2010- 
268246, PIEF-GA-2011-300984, the EU program "Thales" ESF/NSRF 2007-2013, and by 
the European Science Foundation "Holograv" (Holographic methods for strongly coupled 
systems) network. It has also been co-financed by the European Union (European Social 
Fund, ESF) and Greek national funds through the Operational Program "Education and 
Lifelong Learning" of the National Strategic Reference Framework (NSRF) under "Funding 
of proposals that have received a positive evaluation in the 3rd and 4th Call of ERG Grant 
Schemes" . 



- 28 - 



A Integral identities 



We collect the relevant integral identities used in preceding sections. Many of them are 
derived and listed in Appendix of [53]. 

A.l Integrals for the charge density 

In the following, we assume that x £ (a, 6) and T ^ (a, 6). The notation V denotes the 
principal value. We have 



dx 



TT 



V 



a (r±x)V(6-x)(x-a) ^iT±a){T±b) ' 

: 

a {x-y)y^{b-y){y-a) 

Ja [y-x)[y±T) \ T-x 



V 



V{b- y){y - a) 

y -X 

V{b - y){y - a) 

y±T 



dy = + b — 2x) 



dy = -TT 
^ 2 



a + 6 ± 2r =F 2y/{T±a){T±b) 



V 



1 



y-a 



dy = IT 



la {y-x)\l 1-y 
A. 2 Integrals for the free energy 

Let us present the integral identities used in Section 3.3. 

The integrals 

We compute the following integrals: 



^ Jog{T- x)^{b-x){x-a) 



dx 



(l±x) 

log(T — x) (6 — x)(x — a) 
^(6-x)(x-^ 1^^^^ 



^x- 



^/(b- x)(x - a) 



^{x-a) + {l±b) 



^1 dxj"^ — ^log(r-x) dx + (l±6) 



b — x 



(l±a)(l±6) 



log(r - x) 



a (1 ± x) Y^(6^^3?)(x'^^aJ 



(l±a)(l±6) 
1 ±x 

' log(T - x) 
a \/{b- x)(x - a) 

dx 



dx 



^Ai + (1 ± b)A2 - (1 ± a)(l ± b)A^ , 



(A.l) 
(A.2) 
(A.3) 

(A.4) 
(A.5) 
(A.6) 



(A.7) 
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where, using (6.1) and (6.5) of [53], we have 
Ja y/[b - x)[x - a) 

dx 



^3 



['> logT-x) fVT^+Vr^ 

I —j^^^=^^^= dx = 27r log 
a \/(h-x){x-a) 

^ log(T - x) 

a (1 ± x)i/(6 - x){x - a) 

( ±{T ± 1)2 ^ (^(r-a)(r-6) T ^{l±a){l±b)Y ^ 



(A.8) 



= , log , 

V(l±a)(l±6) 

and using (6.1) and (6.4) of [53], we obtain 



(^/rT^+\/rx6)2 



(A.9) 



X — a 

dx\ I log(r — x) dx 

b — X 



x — a 



^ xiog(r-x) 



log(T — x) dx 



dx — a 



log(T — x) 



IT 



^y {b - x){x - a) J a \/{b- x){x - a) 

hb + a)-T+^{T-a){T-b) 



dx 



+ {b- a) log 



(A.IO) 



The integrals 

We compute the following integrals: 



-'2 



dx log(l — x 



Y^(6-x)(x-a) 



l±x 
(A")t=i , 

'dxlog(l + x)^ -"^^"-"^ 

1 ± X 



= {Ii )T=l,a^-a,b^-b 

The integral 

We compute the following integral: 



(A.ll) 



(A.12) 



, log(r-x)A/(6-x)(x-a) 
dx 



J a 

I 



dx 



(T-x) 
log(T — x) (6 — x)(x — a) 
^(b-x)(x-a) (T-x) 



^/{b^^x){x^^ 



(x - a) + (T - 6) + 



{T-a)iT-b) 



dx^ 1^ — - log(r - x) dx + (T - b) 
b — X 



x-T 
^ log(r - x) 
^/Jfi^^xji^r^ 



dx 
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{T-a){T-b) 



log(T — x) 



la {T - x)^y{b - x){x - a) 
=: Ai + iT- b)A2 -{T-a){T- ftjis , 



dx 



(A.13) 



where 



A. 



log(T — x) 



a {T - x)y/{b- x){x - a) 



dx 



2tt 



log 



1 1 

+ 



(A.14) 



The integrals 

We compute the following integrals: 



T± 



Ax log(l lb x) ^ 



T — X 



^/{b-x){x-a) 



{x-a) + {T-b) + 



iT-a){T-b) 



x-T 



X — a 
b — X 



log(l ± j;)dx + (T - 6) 



log(l lb x) 



a \/ {b - x){x - a) 

-(T-a)(T-,) f '°^"^^) 

J a {T - x)yj{b — x) {x — a) 

+ (T - b)Bf -{T-a){T- b)B^^^ , 



dx 



dx 



where, from identities (6.6) and (6.9) of [53], we have 



B 



(2) 
± 



B 



(3) 



y^{b- x)(x - a) 

^ , log(l lb X 

dx 



27r log 



(A.15) 



(A.16) 



, = log , 

^{T-a){T-b) 

and from (A. 10), we obtain 

cb 



(T — x) ^J{Ji^^a^j{x^^^~a) 

( ±{T ± 1)2 ^ (^(T-a)(r-6)T ^(l±a)(l±6))^ 



vr 



(A.17) 



5. 



(1) 
± 



dxA / log(l ± x) 

b — X 



^(fe + a)±lTV(l±a)(l±&) + (fe-a)log( ^^+^^ 



(A.18) 
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B Method of conformal mappings and mirror charges 



In this section, we present an alternative method to find the equihbirium charge density 
cr(x) of the Coulomb gas, which makes use of conformal mappings and mirror charges. This 
method could be used to derive the results discussed in the Sections 3-5. The strength of this 
method is that the free energy, which is required to derive the saddle-point approximation for 
the Hilbert series, can be calculated without solving the equilirium condition explicitly. As 
an example, we apply it to the calculation of Section 5, and solve the free energy at leading 
order. 

B.l Solving the complex potential 

We construct the complex potential 

1 

^(x) = -niog{x - T) - J dy ao{y)logix-y) , (B.l) 
where the LO charge density ao{x) is fixed by requiring that the real part Re <I>(x) is constant 

Re Mx)] = ^Ao (B.2) 

on the support of the charge density [a, 1]. By comparing (B.6) to the definition (5.5) we see 
that (B.2) indeed matches with the LO constant of Section 5. The lower limit a (with 
— 1 < a < 1) is to be determined subsequently. 
We start by applying the conformal map 



X ^ z = fix):= (B.3) 
V X — a 

which maps the exterior of the "conductor" [a, 1] in the compactified complex plane to the 
right half-plane. The point charge —TZ/2 at x = T is mapped to z = Z, where 



Z:=/(T) = ^^--^ and 0<Z<1. (B.4) 

Notice also that there is a point particle with the charge TZ/2 — 1 at x = oo, which is mapped 
to z = 1. The solution for which the real part of the potential is constant on the image of the 
conductor, the imaginary axis, is now obtained by placing mirror charges in the left half-plane 
and reads^'^ 

$(z) = l-niog ('-^) - (In - 1) log ('-^) + $0 . (B.5) 



2 ° \z + Z J \2 J ° \z + l 



^^The complex potential of a point charge e located at the point zq in two dimensions is given by ^{z) 
-elog(z - zo). 
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This solution is the unique one having the desired analytic structure, up to irrelevant branch 
choices and the value of the constant $o- Mapping back to the x-plane, we find the solution 



«>(x) = CD (/(x)) 
1 



x-l 
x—a 



o^log 

Z I x—1 



T-1 
T-a 



+ 



T-1 
T-a 



zTZ-l] lo. 



x-l 



x-l 



+ ^0 . 



+ 1 



(B.6) 



B.2 The free energy 

Evaluating the real part of Eq. (B.6) on the conductor, say at x = 1, we find that the constant 
in (B.2) is 



^0 = 2 Re [$(x = 1)] = 2 Re ^>o . 



(B.7) 



The value of Re $o, and thus also the constant A, is found by comparing the large x asymp- 
totics of the solution (B.6) to the definition (B.l). In this limit, Eq. (B.l) gives 



Re$(x) = ( -7^ - 1 ) log(x) + O 



X — )• CXD 



(B. 



whereas Eq. (B.6) gives 



Re$(x) = -71 log 



T-1 
T-a 



1 + 



T-1 
T-a 



zTl-l] lo, 



-4^ +2^° + ^ 



X — >• oo . 



Equating these two expressions, we have 



Ao = (n- 2) log 



I- a 



71 log 



T-1 
T-a 



1 + 



T-1 

T-a 



(B.9) 



(B.IO) 



On the other hand, at LO the integral / of (2.19) is given by 



In = Tl log(T — x) ao(x) dx = —TZ lim Re 

Ja 

as seen from the definition (B.l). Inserting here the solution (B.6) gives 



^>(x) - -7^1og(x -T) 



(B.U) 



7^ 



(27^ - 2) log 



1 - a 



niog [{T - l)(r - a)] + (2 - 2n) log ■ 



T-1 
T-a 



1 + 



T-1 
T-a 



(B.12) 



The result for the LO free energy is thus proportional to the sum of (B.IO) and (B.12). 
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B.3 The solution for the gapless phase 

The solution for the gapless phase is now given by setting a = —1 in (B.IO) and (B.12): 

Aq = 21og2 - 7^1og(2^) , (B.13) 
/o = -7^1og(2t) +7e2log(l , (B.14) 

where we used (2.3) to replace T by t. Observe that the results coincide with (5.18) and (5.19). 

B.4 The solution for the gapped phase 

In order to determine the solution for the gapped phase, we need to determine the value of a. 
To this end, it is necessary to study the behaviour of the charge density near its lower edge. 

Since conformal mappings conserve charge, it is useful to study the density in the z- 
plane, where x = a is mapped to z = oo and x = 1 is mapped to z = 0. We interpret the 
left half-plane as the conductor with a constant potential, whereas in the right half-plane the 
potential is given by Eq. (B.5). The charge density is the discontinuity of the electric field at 
the surface, 



d$(z) 



dz 

where C = Imz and Z is given by (B.4). 



z[n - (7^ - 2)z] + einz - (7^ - 2)] ^ ^^^^^^ 



2^(l + C2)(C2 + Z2) 



Conditions for gap formation Let us assume that we are in the gapless phase (a = —1). 
The gap is formed when the charge density vanishes at x = — 1, which is mapped to C = ±00. 
As — )■ ±00, the charge density behaves as 

d{C) ^ C'^IT^Z - {n-2)] , C^oo (B.16) 

Then, requiring (B.16) to be non-negative yields 

where we used the definitions of Z and T from (B.4) and (2.3). This equality results in the 
same conditions for the gapless phase as in Sections 4 and 5,7^<2or0<t<tc = l/('^ — 1)- 
In the opposite case, the system is in the gapped phase. 

The lower endpoint of the conductor in the gapped phase To solve the location of 
the lower endpoint a = of the charge density in the gapped phase (with —1 < oq < 1), 
we require that the leading term of the charge density (B.16) vanishes at the endpoint. This 
gives 



from which the value of ag can be solved: 



7^2 - 4(7^ - i)r ,^ , 
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The free energy in the gapped phase The free energy in the gapped phase can now be 
obtained by inserting the solution (B.19) to equations (B.IO) and (B.12). It is easy to check 
that the result agrees with (5.20) and (5.21). 

B.5 Charge densities 

As a byproduct of the method, we can also solve the charge density on the x-plane by applying 
the conformal map (B.3) to (B.15). Let x be a coordinate on support [a, 1] of the density 
a{x). We obtain^^ 

ao{x) := -- lim $'(x + ie) = 2i lim f'{x + ie) a{-if{x + ie)) . (B.20) 

Vr e-5>0+ e-S>0+ 



The explicit results for the charge densities can readily be computed: 
£70 (x) 



X £ [qq, 1] (gapped phase) 
These results are in accordance with (5.9) and (5.16). 
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